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O ■ ABSTRACT 
(N 

Context. M81 is an ideal laboratory to investigate the galactic chemical and dynamical evolution through the study of its young and 
old stellar populations. 

Aims. We analyze the chemical abundances of planetary nebulae and H ii regions in the M81 disk for insight on galactic evolution, 
and compare it with that of other galaxies, including the Milky Way. 

Methods. We acquired Hectospec/MMT spectra of 39 PNe and 20 H ii regions, with 33 spectra viable for temperature and abundance 
analysis. Our PN observations represent the first PN spectra in M81 ever published, while several H ii region spectra have been pub- 
lished before, although without a direct electron temperature determination. We determine elemental abundances of helium, nitrogen, 
oxygen, neon, sulfur, and argon in PNe and H ii regions, and determine their averages and radial gradients. 

Results. The average 0/H ratio of PNe compared to that of the H ii regions indicates a general oxygen enrichment in M81 in the 
^ I last -10 Gyr The PN metallicity gradient in the disk of M81 is Alog(O/H)/ARG=-0.055+0.02. Neon and sulfur in PNe have a radial 

distribution similar to that of oxygen, with similar gradient slopes. If we combine our H ii sample with the one in the literature we find 
a possible mild evolution of the gradient slope, with results consistent with gradient steepening with time. Additional spectroscopy is 
needed to confirm this trend. There are no Type I PNe in our MSI sample, consistently with the observation of only the brightest bins 
Q , of the PNLF, the galaxy metallicity, and the evolution of post-AGB shells. 

» I , Conclusions. Both the young and the old populations of MSI disclose shallow but detectable negative radial metallicity gradient, 

^ ' which could be slightly steeper for the young population, thus not excluding a mild gradients steepening with the time since galaxy 

^ ' formation. During its evolution MSI has been producing oxygen; its total oxygen enrichment exceeds that of other nearby galaxies. 

Key words. Galaxies: abundances, evolution - Galaxies, individual: MSI - Planetary nebulae - H n regions 

> ' 

1. Introduction stars (the so-called o'-elements) probe the environment at the 

epoch of the formation of the PN progenitors. Oxygen, neon, sul- 

O ■ The study of chemical evolution of galaxies hinges on the con- f^^^^ ^^^^^ ^he most part, manufactured in high mass 

; frontation between sets of evolutionai-y chemical models and ai-- ^^^^^ (m>8 M^), thus their concentration in PNe across galactic 

^ rays of observational data. Comparison of models with rehable ^^^^^ ^^^^^ metallicity gradient in remote epochs. The metal- 

O abundance determinations are needed to understand the frame- jj^j^y gradients of PNe, compared to gradients of young stellai- 

O work of galaxy evolution, and to test star formation history, stel- populations (such as those given by H ii regions) allow the study 

T-H . lar nucleosynthesis, and the mechanisms of galaxy assembly and chemical evolution of galaxies. Galactic metallicity gradi- 

^ , formation. ^^^^ measured in terms of Alog (X/H)/ARg dex kpc"', where 

. ^ Galactic metalUcity gradients, and especially their time evo- x/H is the abundance of element X with respect to hydrogen and 

lution, are important probes of galaxy formation and evolution, j^^ jg jj^g galactocentric distance 
H Observationally, most studied galaxies display a radial metallic- 

. 5t , ity gradient with a negative slope, indicative of higher metallic- Planetary nebula gradients in the Milky Way indicate that the 

ity towai-d the inner galaxy. Time variation of gradient slopes metallicity decreases inside out by 0.01 to 0.07 dexkpc"' (e. g., 

in galaxies is hard to come by, given the lack of reliable time- Maciel & Quireza 1999, Henry et al. 2004, Stanghellini et al. 

tagging metallicity probes. From the viewpoint of theory, dif- 2006, Perinotto & Morbidelli 2006). The most recent of these 

ferent models predict opposite time evolution of the metallicity studies seem to agree that the PN metallicity gradients are rather 

gradient, showing the sensitivity of this physical parameter to flat, i.e., their slope is not steeper than —0.02 dex kpc"', and 

the adopted parameterization of the processes (e.g., Magrini et only slightly shallower than that derived from H ii regions, im- 

al. 2007). A key factor to forwarding our knowledge in this field plying that the Galactic metallicity gradient does not vary very 

is the acquisition of high quality spectra to obtain abundances much within the time frame of Galactic evolution. The metal- 

for a variety of probes, in a wide selection of different galaxies, licity gradients based on Galactic PNe are hindered by the un- 

including our own, to cover as much as possible the pai-ameter certainties in PN distances, by the interstellai- absorption in the 

space of initial masses, metallicities, and galaxy type. direction of the Galactic center, and by the paucity of PN ob- 

The planetary nebula (PN) abundances of elements that are servations in the direction of the Galactic anti-center. These fac- 

invariant during the evolution of low- and intermediate-mass tors, together with the interest of studying metallicity gradients 
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of Other spiral galaxies of different morphological types and be- 
longing to different groups, prompt us to investigate the PN (and 
H n region) metallicity gradients in the M8 1 disk. 

M8 1 is the largest member of the nearest interacting group of 
galaxies at a distance of 3.63+0.34 Mpc (Freedman et al. l2001l) . 
Its two brightest companions, M82 and NGC 3077, are located 
within a projected distance of 60 kpc. The disturbed nature of the 
system is evident from the observations in the 21 cm emission 
line of H i, showing extended tidal streams and debris between 
the galaxies (Gottesman & Weliachew'1975', Yun et al. I1994I ). 
The simulation of Yun d 19991 1 explains many of these tidal fea- 
tures as being the result of close encounters between M81 and 
each of its neighbors ~200-300 Myr ago. The most prominent 
debris (e.g., Sabbi et al. 2008, Weisz et al. 120081) are consid- 
ered tidal dwarf galaxies, new stellar systems that formed in gas 
stripped from interacting galaxies. 

The a-element abundances of PNe and H ii regions in M8 1 
have remained elusive, with the exception of 17 H ii regions 
(Garnett & Shields 1987, GS87). With the PNe (Jacoby et al. 
1989, Magrini et al. 2001) and H ii regions (Lin et al. 2003) 
locations, and their de-projected (PA=157'^ and inclination=59", 
Kong et al. 2000) galactocentric distances available we are in the 
position to push forward a determination of abundance gradients 
from PNe and H ii regions in the M8 1 disk. The location of M8 1 
within its group will allow us to study the effect of galaxy in- 
teractions and mergers in the evolution of a galaxy and of its 
abundance gradient. 

In order to sample the M81 disk to the necessary S/N for 
abundance analysis we need spectroscopy with a 6-8m class tele- 
scope. MMT/Hectospec medium-resolution spectroscopy offers 
the unique opportunity to study both the PNe and the H n re- 
gion populations with the same setup, during the same night, and 
adopting the same analysis techniques, thus avoiding most biases 
that commonly affect the comparison of chemical abundances in 
population of different ages (e.g. Fe/H from old RGB stars is dif- 
ficult to compare with O/H from the present time H ii regions). 
We present here the results from the study of a sample of bright 
PNe in M81, with the goal of obtaining the metallicity and ex- 
plore their gradients within the M8 1 disk. Their properties as a 
group are discussed, and compared with those of H ii regions 
from both our MMT observations and from the literature. The 
combined properties of M8 1 PNe and H ii regions are also com- 
pared with their homologs in other nearby galaxies and in the 
Milky Way. 

In §2 we describe the observation and analysis techniques, 
and give the measured fluxes and calculated abundances. §3 in- 
cludes the analysis of the derived abundances and the determi- 
nation of the metallicity gradients, §4 presents the discussion of 
our results, and the conclusions are in §5. 



2. The MMT observations: data reduction and 
analysis 

2.1. Observations 

We obtained the spectra of PNe and H ii regions in M8 1 with the 
MMT Hectospec fiber-fed spectrograph (Fabricant et al. 2005). 
The spectrograph was equipped with an Atmospheric Dispersion 
Corrector and it was used with a single setup: 270 mm ' grat- 
ing at a dispersion of 1.2 A pixel The resulting total spectral 
coverage ranged from approximately 3600 A to 9100 A, thus in- 
cluding the basic emission-lines necessary for the determination 
of physical and chemical properties. 



The PNe were selected among those observed with the INT 
through with [O iii]. Ha, and Stromgren Y filters (Magrini et al. 
2001), which allowed us to cover the whole field of M81, and 
to identify a large number of disk PNe at large galactocentric 
distances. An earlier survey by Jacoby et al. (1989) was limited 
to the M81 bulge. It is worth noting that the INT images allow 
us to have accurate positions (<0.5") and finding charts for the 
selected PNe and H ii regions, a crucial information given the 
crowding of the inner and spiral arm fields of M8 1 . 

The first target selection was performed by scaling the PN 
brightness in [O m] with those of M33 PNe, whose medium- 
resolution spectra have been observed by us using the identi- 
cal technique (Magrini et al., 2009, 2010, hereafter M09, MIO). 
We also endeavored to have PN targets at several galactocen- 
tric distances, in order to characterize the PN metallicity gradi- 
ent. The instrument deploys 300 fibers over a 1 degree diame- 
ter field of view and the fiber diameter is- 1.5" . The projected 
size of M8 1 on the sky is smaller than the whole field of view 
of MTT/Hectospec, and we took advantage of this by placing 
some of the fibers in the outermost periphery of M81, where 
some H ii regions and few PNe, both belonging to the outer disk 
and to the intra-group population, have been identified. In ad- 
dition, since the available catalog of H ii regions by Lin et al. 
(2003) is based on Hff photometry, our INT [O m] observations 
have been useful to select [O iii] -bright H ii regions essential for 
chemical abundance determination. As stated by Magrini et al. 
(2001), we expect the misclassification of H ii regions into PNe 
and vice-versa to be lower than~3% of the whole sample. 

The observations were carried out in queue mode in four runs 
during the months of November and December 2008. Each of 
the observing runs consisted of 3 to 5 exposures of 1800s each, 
for a total of 16 exposures (or 8 hours of observation). We used 
the same fiber setup during all observing runs, and thus PNe and 
HII regions were observed within the same conditions, including 
the same exposure time. Priority was given to PNe when plac- 
ing the fibers, since a set of H n region spectra already exists in 
the literature (GS87). A large number of fibers were devoted to 
sky measurements, to ensure that sky spectra were obtained in 
several positions around each target. 

Cosmic rays were removed with the appropriate hec- 
tospec. hcosmic routine in the hectospec analysis package. 
Cosmic rays are detected via subtraction of multiple images, 
flagging the high or low pixels. The flagged pixels are interpo- 
lated over, and the resultant images are combined by average. 
This method avoids problems with clipping that other programs 
suffer from. In order to perform the sky subtraction we first in- 
spected the sky spectra and eliminated those with unusual sig- 
nals. Then we averaged the six spectra closer (in fiber location) 
to the targets, scaled the result to match the object spectrum, and 
subtracted. 

The spectral calibration was achieved with the standard star 
Hilt600 (Massey et al. 1988), observed during the first night. 
This star was chosen by its excellent spectral coverage compared 
to the spectra of the PNe. It is worth noting that our analysis does 
not hinge on flux calibration, rather spectral calibration and flux 
ratios, thus this approach is sufficient to the science goals. Other 
stars were observed during the run nights, but the excellent spec- 
tral coverage of Hilt600 makes it the ideal choice for this type of 
calibration. Variations in sky conditions are taken into account 
since sky subtraction was performed before flux calibration. 

In Table 1 we give the IDs, equatorial coordinates, and equiv- 
alent [O III] magnitude of the planetary nebulae and H ii regions 
observed with MMT medium-resolution spectroscopy. All PN 
IDs are from Magrini et al. (2001), except for PN4 and PN5, 
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which are identified PNe with unpubHshed positions (Magrini, 
private communication). There are 39 PNe and 20 H ii regions 
whose spectra have been acquired and analyzed. With the excep- 
tion of PN 49m and PN 79m, whose spectral lines are below the 
detection limit, we detect at least the major spectral lines in each 
target. 

In Table 2, available online in its completeness, we give for 
each target the ion and wavelength (columns 1 and 2) of the ob- 
served emission line, then the relative observed flux (column 3), 
its uncertainty (column 4), and the line intensity (column 5) ob- 
tained with the extinction correction given as a header for each 
target. All fluxes and intensities are normalized for Fh/3=100 and 
Ih/3=100. Both PNe and H n regions are listed in Table 2. Target 
names are given for each target at the head of the correspond- 
ing flux list. The emission-line fluxes were measured with the 
package SPLOT of IRAlfl Errors in the fluxes were calculated 
taking into account the statistical error in the measurement of 
the fluxes, as well as systematic errors of the flux calibrations, 
background determination, and sky subtraction. Observations of 
same objects during diff'erent runs were stacked together and the 
resulting spectra have been used for the analysis. 

We derived c(Hy6) from the Ha/Hfi flux ratios, as described 
in M09, and used the Galactic extinction curve by Mathis (1990) 
to correct the observed flux for interstellar extinction. In a few 
PNe the derived intensity of the Hy line is more than 30% off 
the theoretical prediction (14340=0. 466 for Te= 10,000 K, case B, 
Oserbrock & Ferland 2006). The reason for this offset could be 
an uncertain c(Hy6). The average c(Hj6) calculated for our PNe 
is 0.7+0.3, a reasonable value if we consider the average c(Hj6) 
of Galactic PNe (from Cahn et al. 1992, <c(H;6)>=0.95) and the 
relative metallicities of the Galaxy to M8 1 . We will recall these 
outliers later in the discussion and indicate them with different 
symbols in the Figures. 

To our knowledge, this is the first attempt in the literature to 
detect PN emission line fluxes in M8 1 thus we cannot extend our 
discussion to comparison with previous work. The spectral qual- 
ity is excellent. In Figure 1 we show a few examples of our M8 1 
PN spectra, where the plasma and abundance diagnostic lines are 
well evident, and the S/N obtained with our observations makes 
us confident on the data quality. 

Several H 11 regions have been already observed by GS87, 
where the Authors refer to the catalog by Hodge & Kennicutt 
(1983) for coordinates. Since the latter are given in image co- 
ordinates rather than Galactic coordinates, it has been hard to 
eliminate from our sample of H 11 regions those that have been 
already analyzed by GS87 or comparing targets from the two 
sets. But by comparing the position visually, we did not find any 
obvious overlap between the two samples, thus cannot compare 
the results. 

2.2. Plasma diagnostics 

We perform plasma diagnostics by using the 5-level atom model 
included in the nebular analysis package in IRAF/STSDAS 
(Shaw & Dufour 1994), consistently with our work in M33 
(M09, MIO). To calculate densities we use the doublet of 
the sulfur lines [S ii]/l/l6716,6731, while for the electron 
temperatures we use flie ratios [O iii]/i4363/(/i5007 -H/i4959), 
[N ii]^5755/(/l6548 + A6584), and [S iii]/l63 12/(^9059 -1- 



' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation 



Table 1. Observing log 



Id. 




CJ2000 


6nooo 




m,|5007 


(1) 




(2) 




(3) 






(4) 




PN 


3m 


9:54 


26.54 


69 


10 


30.5 


24 


41 


PN 


4 


9:54 


47.50 


69 


35 


42.0 


24 


15 


PN 


4m 


9:54 


27.66 


69 


13 


55.0 


24 


45 


PN 


5 


9:54 


53.80 


69 


08 


01.0 


23 


87 


PN 


6m 


9:54 


36.70 


69 


17 


00.5 


24 


39 


PN 


7m 


9:54 


38.52 


69 


10 


36.8 


24 


53 


PN 


8m 


9:54 


45.50 


69 


08 


06.6 


24 


06 


PN 


9m 


9:54 


46.52 


69 


05 


39.0 


23 


66 


PN 


11m 


9:54 


53.15 


69 


13 


32.2 


24 


51 


PN 


14m 


9:54 


09.28 


69 


11 


31.0 


24 


75 


PN 


15m 


9:55 


01.01 


69 


12 


04.3 


24 


45 


PN 


17m 


9:55 


10.09 


69 


02 


37.1 


24 


24 


PN 


28m 


9:55 


13.62 


69 


07 


10.3 


23 


91 


PN 


29m 


9:55 


14.78 


68 


55 


30.0 


23 


56 


PN 


32m 


9:55 


15.70 


69 


08 


39.8 


24 


23 


PN 


33m 


9:55 


15.96 


69 


13 


31.7 


24 


30 


PN 


38m 


9:55 


19.47 


68 


54 


45.7 


24 


39 


PN 


41m 


9:55 


02.58 


69 


06 


26.3 


23 


77 


PN 


45m 


9:55 


21.81 


69 


07 


50.6 


23 


86 


PN 


49m 


9:55 


23.38 


69 


02 


50.9 


23 


66 


PN 


63m 


9:55 


28.25 


69 


01 


26.5 


24 


26 


PN 


70m 


9:55 


30.62 


69 


07 


56.3 


23 


66 


PN 


79m 


9:55 


34.58 


69 


09 


00.9 


24 


49 


PN 


93m 


9:55 


39.07 


69 


05 


42.1 


23 


48 


PN 


117m 


9:55 


50.11 


68 


59 


48.4 


24 


04 


PN 


121m 


9:55 


51.46 


69 


05 


14.9 


23 


56 


PN 


127m 


9:55 


52.93 


69 


10 


52.1 


24 


60 


PN 


128m 


9:55 


53.02 


69 


06 


15.2 


24 


19 


PN 


136m 


9:55 


06.11 


68 


57 


18.5 


24 


60 


PN 


137m 


9:55 


07.29 


69 


12 


10.4 


24 


47 


PN 


147m 


9:56 


20.74 


68 


58 


26.8 


24 


41 


PN 


149m 


9:56 


27.09 


69 


05 


26.4 


23 


66 


PN 


153m 


9:56 


28.70 


69 


00 


32.8 


23 


56 


PN 


157m 


9:56 


31.99 


69 


03 


31.3 


24 


09 


PN 


158m 


9:56 


32.84 


69 


03 


58.8 


24 


23 


PN 


160m 


9:56 


04.48 


69 


06 


36.0 


24 


31 


PN 


166m 


9:56 


08.41 


69 


08 


37.8 


24 


49 


PN 


167m 


9:56 


08.54 


69 


03 


56.0 


24 


53 


PN 


169m 


9:56 


09.87 


68 


59 


56.0 


24 


39 


Hii 


4 


9:54 


34.67 


69 


05 


52.2 






Hii 


5 


9:54 


36.1 


69 


06 


30.5 






Hii 


31 


9:54 


43.62 


69 


03 


29.9 






Hii 


42 


9:54 


47.66 


69 


02 


18.7 






Hii 


72 


9:54 


52.93 


69 


08 


51.8 






Hii 


78 


9:54 


53.84 


69 


01 


08.8 






Hii 


79 


9:54 


54.28 


69 


10 


56.4 






Hii 


81 


9:54 


54.53 


69 


09 


30.5 






Hii 


121 


9:55 


02.87 


69 


04 


49.3 






Hii 


123 


9:55 


03.64 


69 


10 


54.6 






Hii 


133 


9:55 


04.98 


69 


09 


53.6 






Hii 


201 


9:55 


22.82 


69 


09 


52.9 






Hii 


213 


9:55 


26.01 


69 


11 


58.9 






Hii 


228 


9:55 


32.66 


69 


11 


55.4 






Hii 


233 


9:55 


34.59 


69 


08 


31.1 






Hii 


249 


9:55 


37.61 


69 


12 


03.6 






Hii 


262 


9:55 


41.42 


69 


11 


12.2 






Hii 


282 


9:55 


46.50 


68 


59 


24.6 






Hii 


325 


9:55 


58.48 


69 


09 


9.9 






Hii 


328 


9:55 


59.29 


68 


59 


37.6 






Hii 


348 


9:56 


03.23 


69 


07 


38.3 






Hii 


352 


9:56 


03.79 


69 


09 


37.8 






Hii 


384 


9:56 


08.33 


69 


05 


14.8 






Hii 


403 


9:56 


11.78 


69 


07 


44 
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Fig. 1. Sample of medium-resolution spectea of M81 PNe. 



/19532). In some cases, even if the sulfur lines are available, 
their ratio is outside the diagnostic ranges for density calcula- 
tion (e.g., Stanghellini & Kaler 1989), and we assumed a den- 
sity of lO^cm"^. We calculate the medium- and low-excitation 
temperatures respectively from the [O m] and [S m], and from 
the [N ii] line ratios, respectively (see also Osterbrock & Ferland 
2006,§5.2). 



2.3. Abundance analysis 

Of the 39 PNe with emission line spectra, [O m] and [N n] 
temperature diagnostics are available respectively for 10 and 12 
PNe, with only 3 PNe in common. In the analysis of the M33 
PNe (M09) we were able to infer the electron temperature for 
targets whose plasma diagnostics were not available by using 
the correlation between the electron temperature and the He n 



line intensity. In the case of M81 none of the PNe show the 
He II emission; as a consequence, abundance analysis was pos- 
sible only for the 19 PNe where electron temperature was diag- 
nosed directly. There are also 14 H n regions whose temperature 
analysis was based on auroral lines, and their plasma and abun- 
dance analysis was performed in a similar way than for the PNe. 
Electron temperature diagnostics for more than one excitation 
state was available for 5 of the 14 regions. 

The plasma diagnostics are given in Table 3, available online, 
together with the ionic abundances of the PNe and H n regions. 
Column (1) gives the diagnostics, and column (2) the value as 
determined from our analysis. The targets are listed sequentially, 
as in Table 1, and the first hue for each target gives the target 
name. 

The abundance analysis of the PNe and the H ii regions 
was performed following the prescription given by M09, MIO; 



Stanghellini, L. et al.: M81 draft 



Table 2. Observed Fluxes and Line Intensities 



Ion 


^(A) 


Fa 






(1) 


(2) 


(3) 


(4) 


(5) 


PN 3m, C;5=0.502+0.027 










[On] 


3727 


373.9 


58.6 


517.3 


Hy 


4340 


50.0 


3.5 


58.0 


W 


4861 


100.0 


4.5 


100.0 


[Ora] 


4959 


35.4 


3.8 


34.4 


LOraJ 


DUU/ 


104./ 


J.Z 


99.9 


[Nil] 


5755 


1.3 


0.3 


1.1 


He I 


5876 


9.1 


0.3 


7.2 


[Nn] 


6548 


59.8 


3.5 


42.2 


Ha 


6563 


412.0 


6.9 


290.1 


[Nn] 


6584 


142.5 


4.5 


100.0 


[Sii] 


6717 


60.3 


3.8 


41.4 


[Sii] 


6731 


44.3 


3.4 


30.4 



Table 3. Plasma Analysis and Abundances 



(1) 


(2) 


Name 


PN3m 


r,[Nii] 


8893 


Ne[S II] 


100 


Hei/H 


0.051 


He/H 


0.051 


[On]/H 


3.320x10-" 


[O m]/H 


5.440x10-^ 


ICF(O) 


1.0 


0/H 


3.864x10-'' 


12+log(0/H) 


8.587 


[N n]/H 


2.916x10-^ 


ICF(N) 


1.164 


N/H 


3.394x10-"^ 


12+log(N/H) 


7.531 


[S n]/H 


2.250x10-* 


ICF(S) 


1.001 


S/H 


7.066x10-* 


12+log(S/H) 


6.849 



abundances of all elements except helium where calculated with 
the ionization correction factors (ICFs) given in Kingsburgh 
& Barlow (1994) for the case were only optical lines are de- 
tected. Where both Te[0 III] and Te[N II] were available we 
used Te[0 III] to derive the abundances of medium and high 
excitation ions, and Te[N II] for the low excitation ones. In sev- 
eral H II regions the temperature from the auroral line of [SIII] 
were also available, corresponding to transition of intermediate 
excitation. If only one of the electron temperature diagnostics 
was available we use that temperature to derive abundances for 
all ions. Helium abundances were calculated with the formu- 
lation of Benjamin et al. (1999) in two density regimes, and 
taking into account Clegg's (1987) collisional populations, as 
in M09. The sulfur abundances for our PNe are mostly based 
on [S ii] lines, thus highly unreliable. We used the interpola- 
tion to the Galactic PN data by Kingsburgh & Barlow to infer 
[S iii]/[S ii] from the [O iii]/[0 ii] ratios, where possible, and 
obtained an estimate of [S iii], which is observed only in three 
PNe. In order to check whether the sulfur correction is sensi- 
ble, in Figure 2 we plot log(S2+/S+) versus log(02+/0+) for 
the Kingsburgh & Barlow sample. On the same plot we locate 
our M81 PNe with measured [S iii]. The solid line, represent- 
ing S2+/S+ = 4.677 x log(02+/0+)"•'^3^ has the correct slope for 
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Fig. 2. The log(S2+/log(S+) versus log(02+/log(0+) correlation 
(solid line), based on the Galactic PNe from Kingsburgh & 
Barlow (1994, open circles). Our targets are shown with filled 
circles. 



both the Galactic and M81 data, but it overcorrect slightly for 
[S III]. Since the [S iii] lines we observe in M81 are quite faint, 
the relation by Kingsburgh & Barlow (1994) seems reasonable. 

The formal errors in the ionic and total abundances were 
computed taking into account the uncertainties in (i) the ob- 
served fluxes, (ii) the electron temperatures and densities, and 
(iii) c(Hy6). We performed complete formal error propagation for 
a subsample of target, representing a spread of flux error and 
magnitudes, and realized that the errors in the electron tempera- 
ture completely dominate the errors in the ionic abundances for 
all elements except helium. In particular, the flux uncertainty of 
the weaker of the diagnostic lines in the electron temperatures, 
^5755 for Te[N II] and /14363 for Te[0 III], determine the er- 
rors in the abundances. In Table 4 we give the error analysis for 
the PN subsample, including the ones with the minimum and 
maximum flux errors; column (1) gives the PN name, column 
(2) gives the relative flux error for the weaker line in the tem- 
perature diagnostics (in PN 9m we could calculate both Te[N 
II] and Te[0 III], thus we give the two errors), columns (3) and 
(4) give the relative errors in Te[N II] and Te[0 III] respectively, 
and finally columns (5) through (7) give the abundance errors, in 
dex. We derived the errors for all other targets by interpolation 
of the AF^-AA(X) curves. In the case of the helium abundances 
we run the error analysis by using the formulae by Benjamin 
et al. (1999), and found that the factor determining the uncer- 
tainties is the relative error on the helium flux lines. The abun- 
dances, their uncertainties, together with the galactocentric dis- 
tances (see §3.1) are given in Table 5. 

It is worth noting that another source of error that we could 
not estimate explicitly is intrinsic to the ICF method. The un- 
certainties due to the ICF are very moderate for oxygen, but can 
be quite large for the other atoms, especially sulfur, thus these 
abundance uncertainties should be considered as formal errors, 
but might underestimate the real uncertainties. Another source 
of abundance uncertainties is the use of a single electron tem- 
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Fig. 3. log (N/O) vs. He/H; broken Unes show He/H=.125 and 
log(N/O)=-0.3. Type I PNe would be located in the upper-right 
quadrant. 



perature for both low- and high-ionization ions. This regards the 
16 PNe and 9 H n regions where only one choice of the electron 
temperature was available. We have analyzed the sample of tar- 
gets where two electron temperatures were available, and found 
that the average difference between electron temperatures of the 
order of 20%. We tested that the resulting atomic abundance un- 
certainties are lower than those due to flux errors on the auroral 
Unes. 

The sets where we have oxygen and neon in common is lim- 
ited to 6 PNe and 7 H ii regions. There is a clear correlation 
between these two a-elements in both sets, with correlation co- 
efficients Rxy=0.7 and 0.9 respectively for PNe and H n regions. 
This high correlation is expected, since both elements are pro- 
duced in the Type 11 SN environment. The large uncertainties in 
the sulfur abundances make the comparison between oxygen and 
sulfur less meaningful, although the two sets have higher overlap 
than the previously described pair (16 PNe and 14 H n regions in 
common), the correlation coefficient is only 0.4 and 0.3 for the 
two sets. 

Figure 3 shows the log(N/0) vs. helium plot, to classify 
Type I PNe. These are the nitrogen and helium-enriched PNe, 
whose progenitors have likely undergo the third dredge-up and 
hot bottom-burning, and thus are likely to have higher progenitor 
masses (Peimbert & Torres-Peimbert 1983, Marigo 2001). The 
top right portion of the plot, where He/H>0.125 and log(N/0)>- 
0.3, is the locus of Type 1 PNe as determined for the Milky Way 
by Perinotto et al. (2004). Since the metallicity of M81 is sim- 
ilar to that of the Milky Way, we can use the same diagnostic 
plot for M81 PNe. Clearly, there are no Type I PNe among those 
observed by us in M81. In the next section we discuss the popu- 
lation of the observed PNe and its implications for gradients and 
evolutionary studies, and show that indeed we do not expect an 
observed population of Type 1 PNe among our targets. 

In Table 5, at the bottom of each set of targets, we give the 
average abundances, relative to hydrogen (linear for helium, and 
in logarithmic format for the other elements. Uncertainties rep- 
resent ranges of abundances). The average oxygen abundance 



of M81 PNe (2.83+2.01x10"'') is lower that what has been in- 
ferred for the Galactic disk PNe (4.49±2.28xl0-'' for Type II 
PNe, Stanghellini & Haywood 2010) and about 1.6 higher than 
the average oxygen abundance of Type 11 1 PNe in M33 (M09). 

The H II regions observed by us span a small fraction of the 
galaxy radius (5.8 to 9.8 kpc from galaxy center). The average 
O/H value for these regions, (3.42±1.78)xl0"^ should be com- 
pared to the oxygen abundance of PNe in the same radial range, 
which is (1.96+1.27)xl0"'', showing a 0.25 dex oxygen enrich- 
ment in this area of the galaxy. The comparison between PNe 
and H ii region average abundances of other elements is very 
limited, since the sets have 2 or 3 objects in coimnon in the 5.8- 
9.8 kpc region. 

2.4. The population of M81 PNe. 

The PNe whose medium-resolution spectra are presented here 
are selected form the bright end of the planetary nebula luminos- 
ity function (PNLF) of the M8 1 disk. They have two peculiari- 
ties, as a group: first, we do not detect the He n /14686 emission 
Unes in any of them; this is a signature of moderately low cen- 
tral star temperature (less than ~ 100,000 K). Second, none of the 
PNe whose abundance analysis was performed (i.e., whose elec- 
tron temperature has been measured from the emission Unes) is 
of Type I, thus they must be the progeny of stars with tumofF 
masses (Mto) in the~l-2 Mq range. 

The properties of the M8 1 PN sample are important to de- 
fine their nature when we compare the population of M81 PNe 
to those of other galaxies. M81 has similar metalUcity to that 
of the Milky Way (Davidge 2009). In M81, as in most external 
galaxies, the observed PN population belong to the brightest few 
magnitude bins. The nature of the brightest PNe varies notably 
with the metalUcity of the host galaxy. In their runaway evolution 
from the AGB, the PN central stars (CS) reach hot temperature 
while keeping very bright; later, both brightness and tempera- 
ture decline. A PN is visible only when the CS temperature is 
high enough to ionize hydrogen, and the nebula is optically thin 
to the Ha radiation (e.g., StanghelUni & Renzini 2000). At rela- 
tively high metallicities such as those of M81 or the Galaxy, the 
brightest PNe are not those with the hottest CSs (e.g., Villaver 
et al. 2002). In fact, the post-AGB shells with very hot CSs are 
statistically still enshrouded in dust at early phases or their evo- 
lution, thus stiU thick to the Ha radiation. A lower metallicity 
implies a shorter thinning time (Kaufl etal. 1993). Plenty of Type 
I PNe have been observed within the bright end of the PNLF of 
low metallicity galaxies such as M33, the LMC, and the SMC; 
this is not the case in M81, where the high metaUicity prevents 
massive progenitor PNe (i.e., the Type I) to populate the bright 
end of the PNLF. 

To illustrate this point, in Figure 4 we show the luminosity 
function (PNLF) of Galactic PNe in the ^S-equivalent absolute 
magnitude. The Hj8 fluxes are from Cahn et al. (1992), and the 
PN distances from Stanghellini et al. 2008 (see also Stanghellini 
& Haywood 2010). The solid histogram represents the homoge- 
neous Galactic sample of 331 PNe where both distance and HjS 
flux are known, and where spectroscopy has been performed and 
a lower limit or a detection of the He n /14686 line has been ac- 
quired. The shaded histogram is the subsample of Galactic PNe 
with 14686/I/s < 5%. It is clear that at least 80-90% of the bright 
Galactic PNe do not show the He n Unes, in qualitative agree- 
ment with not finding this emission Une in the M8 1 PN sample. 

In order to quantify the expectation of detecting the He n 
emission line in a given PN population of a certain metallic- 
ity we use Kaufl et al.'s (1993) models to derive PN opacity. 



Stanghellini, L. et al.: M81 draft 



7 




40 



30 



20 



10 







Fig. 4. The Galactic PNLF in the light of H/3 for a homogeneous 
sample of PN (unshaded histogram) and for a subsample of PNe 
with no He ii emission lines (shaded histogram, see text). Data 
are from Cahn et al. 1992, Stanghellini & Haywood 2010. 



where the non-evolving term of the opacity,Ti(0), scales di- 
rectly with the dust-to-gas ratio (see their Eq. 10). We also ex- 
press this ratio as expected at different metallicities, Md/Mg ~ 
0.01(O/H)/(O/H)mw (Draine 2009), where the right hand side 
term refers to the ratio of the oxygen vs. hydrogen abundance 
in a given galaxy compared to that of the Milky Way. To obtain 
the ratio for M81 scaled to that of the Galaxy we use the H ii re- 
gion oxygen abundances in M81 (4.99x10"'*) and in the Galaxy 
(5.25x10-'*, Deharveng et al. 2000) to find Md/Mg ~ 0.01, thus 
for a given type of M8 1 PN the thinning time would be similar 
to that of the Galactic homologous. We conclude that the lack of 
He ii-emitting PNe in MS 1 is totally expected, given the Galactic 
PNLF for both He ii-emitting and non-emitting PNe (Figure 4). 

It is interesting to compare the PN population of M81 with 
that of M33 as well. Magrini et al. (2010) obtained an average 
oxygen abundance of <O/H>=2.04xl0 "^for H ii regions in the 
M33 disk, corresponding to a dust-to-mass ratio of~4xl0-"^, or, 
less than half that of M81. By running the Kaufl et al.'s mod- 
els for M33 we obtain (Balmer line) thinning times~1.5 times 
shorter than in MS 1 or the Galaxy, which can account for the rel- 
atively higher frequency of He ii-emitting, bright PNe that where 
observed in M33 by M09. 

We know that the brightest PNe in a given population ex- 
cludes the progeny of the massive AGE stars, i.e., those with 
Mto > 2 Mq (Stanghellini & Renzini 2000). The sample of M81 
PNe under study here thus does not include the progeny of the 
more massive AGB stars. This agrees very well with the obser- 
vational fact that we do not observe any Type I PNe in our M8 1 
sample. 



Fig. 5. Oxygen abundances vs. galactocentric distances for M81 
targets. Top panel: the observed sample of PNe, with formal 
error bars and least square fit represented by solid line. Open 
symbols represent those PNe whose extinction constant might 
be overestimated. Middle panel: The H n region sample from 
MMT spectra (open squares) superimposed to the PN popula- 
tion of the top panel. Bottom panel: The H ii regions from GS87 
(filled squares) superimposed to the MMT H ii regions of the 
middle panel (open squares). The line represents the least square 
fit for the combined sample, as given in Table 6. 



3. Radial metallicity gradients in l\181 

3. 1 . The planetary nebula gradients 

The metallicity gradient of a disk galaxy is generally intended 
as the metallicity variation between the galactic center and the 
periphery, measured radially. These gradients are essential con- 
straints for the chemical evolution models, and, depending on the 
probe, can constrain the galactic composition at different times 
in evolution. Planetary nebulae with low-mass progenitors, such 
as non-Type I PNe, probe the early stages of disk metallicity. 
We could derive the metallicity gradients in M8 1 from oxygen 
abundances, thus probing the a-element evolution through the 
galactic disk. We estimate the distances of our target PNe from 
the center of M8 1 by using a rotation angle of 157° and an incli- 
nation of 59". The error in the target coordinates is very low, and 
the major source of error in the distance calculation is the uncer- 
tainty in the galaxy inclination, which we take as in Connolly et 
al. (1972), +2°, and we propagate to estimate the galactocentric 
distance uncertainties. 

In Table 6, columns (3) and (4), we give the metallicity gra- 
dients calculated for oxygen, neon, and sulfur. Slopes are in dex 
kpc while intercepts are in dex, with the usual notation of 
A(X)=log(X/H)+12. In Figure 5 (top panel) we show the oxy- 
gen gradient from the M81 PNe, where the line represent the 
linear fit. Only PNe are plotted in this top panel, with filled sym- 
bols representing the general PN population, and the open sym- 
bols indicating those PNe whose extinction correction from the 
Ha/Hfi ratio does not well reproduce the line intensity expected 
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Table 6. Metallicity Gradients 
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NpNc 


Alog(X/H)/ARG 


A(X)o 


(1) 


(2) 


(3) 


(4) 


o 


19 


-0.055+0.02 


8.72+0.18 


Ne 


6 


-0.047+0.03 


8.18+0.27 


S 


16 


-0.069±0.01 


7.27+0.10 


X 


Nh ii reg. 


Alog(X/H)/ARG 


A(X)o 





31 


-0.093+0.02 


9.37+0.24 



for the Hy emission line. These open data points might thus have 
larger uncertainties than inferred from the formal error analysis. 

We perform the fits with the routine fitexy (see Numerical 
Recipes, Press et al. 1992) and obtain Alog(0/H)/ARG=- 
0.055±0.02 dex kpc Should we exclude the open symbols in 
the fit we would obtain a very similar result (Alog(0/H)/ARG=- 
0.053±0.02 dex kpc-'). 

In Figure 6 we show the metallicity gradients for sulfur 
(top) and neon (bottom), where the PNe are plotted with filled 
symbols. The solid lines represent the fits obtained with the^- 
texy routines that consider the formal errors both in abundances 
and distances (see Table 6). Fits obtained with the least square 
method are very similar. 

3.2. H II regions in the A(X)-Rq plots 

The bottom panel of Figure 5 shows the H ii region analysis of 
oxygen abundances from this paper (open symbols) and from 
GS87 (filled symbols). Formal errorbars have been applied to 
the data obtained in this paper, while in the case of the GS87 
sample we assume a 0.05% relative uncertainty in the distances, 
and a 0.25 dex uncertainty in the oxygen abundances. It is worth 
noting that we have recalculated the galactocentric distances for 
the GS87 targets by assuming the same distance to M81 than 
we use thorough this paper (d=3.63 Mpc) for homogeneity, al- 
though this does not affect the resulting gradient significantly. 
The composite data sample has been fitted including the distance 
and abundances uncertainties with the fitexy routine, to give a 
gradient (Alog(O/H)/ARG=-0.093±0.02 dex kpc"') not too dis- 
similar, yet steeper, from that of PNe. A direct least square fit 
would give a slightly milder gradient (Alog(O/H)/ARG=-0.07 
dex kpc '). It is worth noting that both the average O/H abun- 
dance in H II regions and the intercept of the linear metallic- 
ity gradient are consistently higher that the respective values for 
PNe, as shown in Table 6, with PNe under abundant with respect 
to H II regions by almost a factor of 2, on average. This is well 
seen in the middle panel of Figure 5, where PNe and H ii regions 
from MMT observations are plotted together 

In Figure 6 we show the sulfur (top) and neon (bot- 
tom) abundances in H ii regions versus their galactocentric 
distances (open symbols), plotted together with the PN data 
(filled symbols). Since the H ii regions observed with the 
MMT span a small fraction of the galaxy radius, metallic- 
ity gradients derived form these species would not be very 
meaningful. Within the distance ranges populated by our 
H II region data we find< Ne/H >h ii / < Ne/H >pn=1.1, and 
< S/H >H ii / < S/H >PN=0.9. Any further comparison is not 
sensible given the limited data at hand, and we prefer to postpone 
a complete discussion of the sulfur and neon abundances in PNe 
as compared to the H ii regions of M8 1 when more H ii region 
spectral data with diagnostic lines relative to these atoms will be 
available. 
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Fig. 6. Sulfur (top panel) and neon (bottom panel) abundances 
from MMT observations vs. galactocentric distances for M81 
PNe (filled circles) and H n regions (open squares). The lines 
represent the hnear fits to the PNe as discussed in the text. 



4. Discussion 

4.1. Galaxy evolution through oxygen abundances 

The O/H abundances are among the most reliable in emission- 
line targets when only optical spectra are available. In the case 
of M8 1 , PNe and H ii regions provide a meaningful compari- 
son between young and old stellar objects, since the techniques 
of data acquisition, analysis, and abundance calculation are the 
same for both classes of objects. In particular, for several H n re- 
gions presented here, those observed with the MMT, the actual 
data acquisition is identical that that of the PNe, thus the com- 
parison among data set is even more compelling. 

It is worth emphasizing that the H n region gradient is not 
based exclusively on our MMT data, which cover a limited radial 
range, but also on the results by GS87. The spectroscopic obser- 
vations of GS87 did not allow a direct measurement of the elec- 
tron temperature. They derived Tg by applying an indirect cali- 
bration made through the photoionization models by Pagel et al. 
(119791) . The errors on the electron temperature are estimated of 
the order of +1200 K, mainly due to the scatter of the measure- 
ments on which Pagel et al. (1 19791 1 based their calibrations. We 
have estimated, using the IRAF ionic task, that an uncertainty 
of + 1200 K in translates in ±0.2 dex in the ionic abundances. 
Thus, the chemical abundances of GS87 result less accurate than 
ours, and may introduce biases in the determination of the gra- 
dient. We conservatively have assumed uncertainties of 0.25 dex 
for the oxygen abundances from GS87. 

As shown in Figure 5, the radial oxygen gradient slope of 
H II regions is mildly steeper than that of PNe, with differences 
compatible with the slope uncertainties. Since there are no Type I 
PNe in our sample it is sensible to assume that the PN progenitor 
(turnoff) masses are smaller than ~2 Mq, implying a population 
older than t ~1 Gyr (Maraston 1998). The lower limit for PN 
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Fig. 7. The slope of the oxygen gradients, in dex kpc ver- 
sus the age of the stellar population considered, for the Galaxy 
(filled circles, adapted from Stanghellini & Haywood 2010), 
M33 (starred symbols, adapted from M09 and MIO) and M81 
(squares, this paper). The zero-age populations represent the 
slopes of the gradients estimated though H ii region oxygen 
abundances, while all other points are slopes of PN gradients, 
see text. Error-bars in the x direction represent the approximate 
age span of the populations, and in the y direction are the esti- 
mated metallicity uncertainties. The zero-age data point for M8 1 
is uncertain, as it is based on the combination of of the H ii re- 
gions observed with the MMT and the GS87 sample 



progenitor mass is around 1 Mq, corresponding to a galactic age 
oft -10 Gyr. 

In Figure 7 we plot the slope of the oxygen radial gradients 
in M81 (squares) for PNe and H ii regions, against the mean 
age of the population considered. H ii regions are placed at t=0, 
while PN ages are inferred from the presumed mass of progeni- 
tors (see Stanghellini & Haywood 2010 for a detailed discussion 
of the different PN types). The horizontal error bar for PNe in- 
dicates the age range of the population, while all vertical error 
bars indicate the uncertainties in the slopes. It appears that the 
slope of the M81 metallicity gradient did not vary very much 
during the past 10 Gyr. A flattening with time, as predicted by 
many chemical evolution models, is not supported by our obser- 
vations. On the contrary, a constant slope, or a mild steepening 
of gradient slopes with time is inferred. 

While uncertainties in the actual slopes are high, it is worth 
comparing the results of M81 with what found in M33 PNe and 
H n regions (M09 and MIO, starred symbols), and the MiUcy 
Way (PNe from Stanghellini & Haywood 2010, H ii regions 
from Deharveng et al. 2000, filled circles). The metallicity gra- 
dients are generally steeper in M8 1 than in either the MiUcy Way 
and M33, but in all three galaxies we observe mild flattening 
with the age of the population, going from steeper for H ii re- 



gions to flatter for non-Type I PNe. It is worth noting that in the 
Galaxy the PN types I, II, and III were independently selected 
to calculate the oxygen radial gradients, while in M33 we were 
able to distinguish between Type I and non-Type I. 

Since the gradient slopes seem to flatten with the age of the 
population considered, there is an indication that the gradients 
themselves are steepening with the galaxy age in all the three 
galaxies considered. It is a mild steepening, and with large un- 
certainties, but it is consistently found in all three galaxies, so it 
is worth exploring it. 

Let us recall that the slope of the H n gradient in M33 ex- 
cluded the central kpc, as explained by MIO. If we were includ- 
ing H II regions in the central kpc we would obtain nearly no 
gradient evolution between H ii regions and Type I PNe in M33. 
There is thus some evolution of the oxygen gradients with time, 
where the slopes vary by 1 x 1 0"^ to~ 3x1 0"^ dex per radial kpc 
in -10 Gyr Stangheflini & Haywood (2010) offer the hypoth- 
esis that portions of the Milky Way at different radial distances 
might have evolved at different rates. The model by Schonrich 
& Binney (2009) illustrates this type of evolution, which could 
hold as well for the other spiral galaxies and for M8 1 . 

In Figure 8 we show the global oxygen evolution in differ- 
ent galaxies, where average oxygen abundances in H n regions 
are plotted against the same quantities in non-Type I PNe. In 
the case of the disk galaxies, both populations cover a reason- 
able radial range and are representative of the whole disks. The 
SMC and LMC oxygen abundance are from Leisy & Dennefeld 
(2006), and the H ii regions in the Magellanic Clouds are from 
Dennefeld (1989). The samples of PNe and H n regions con- 
sidered in the Galaxy and M33 are as described earlier in this 
section. M8 1 PN and H ii regions abundances are from this pa- 
per, both selected consistently within the galactocentric distance 
domain where both sets of targets coexist. The error-bars in the 
lower right of the plot represent conservative data ranges for 
the galaxies considered here. The solar point is from Asplund 
et al. (2005). In this plot M81 stands out as the galaxy with most 
chemical enrichment in the group considered. M33, as found by 
MIO, shows almost no enrichment, and the other galaxies show 
either low enrichment (the Galaxy and the LMC) or slight deple- 
tion of oxygen in the SMC, where the cause could be extremely 
efficient ON cycle at very low metallicity (Karakas & Lattanzio 
2007). The average enrichment of M81 is 0.25 dex. 

The M81 galaxy disk shows indication of chemical en- 
richment through oxygen abundances. It also shows generally 
steeper metallicity gradients than both M33 and the Galaxy. The 
non-type I PNe gradients in M81 is - 0.055 dex kpc while ho- 
mologous slopes are - 0.03 and - 0.02 dex kpc ' respectively in 
M33 and the Galaxy, indicating that whatever process produces 
the metallicity gradients in spiral galaxies, it is more efficient in 
M8 1 and the overall galaxy metallicity is not a factor Similarly 
for H II regions, the M81 gradient slope is - 0.09 dex kpc 
to be compared with - 0.044 and -0.04 dex kpc"' for M33 and 
the Milky Way respectively. Indication of steeper gradients, and 
sizable chemical evolution, are thus the defining characteristics 
of M81 compared with other spirals. The oxygen enrichment in- 
dicates that this galaxy has suffered outflow to a lesser extent 
than the comparison galaxies, and steeper oxygen gradients are 
compatible with this explanation. 

4.2. The chemical evolution of M81: comparison with the 
Milky Way 

It is evident from the above analysis that the metallicity gradient 
slopes in the galaxies examined do not depend on the average 
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Fig. 8. Difference between average oxygen abundance of H ii re- 
gions and non-type I PNe vs. zero-age average oxygen abun- 
dance for M81 (square, this paper), compared to that of M33 
(starred symbol, M09 and MIO), the Magellanic Clouds (LMC: 
triangle, SMC: diamond; Gai-nett 1999, M09), and the Milky 
Way (filled circle, Stanghellini & Haywood 2010, Peimbert 
1999). The solar value represented is from Asplund et al. (2005). 
The M81 datum is derived from PN and H ii regions in the 5.8- 
9.8 kpc range, and does not include the GS87 data.The M81 
galaxy shows considerable metal enrichment. 

Table 7. Properties of M8 1 and the Galaxy 



M81 MW 



Morphological type 




Sa 


Sbc 


Total mass [Mq] 


2.3x10" 


(a) 


1.8-3.7x10" (b) 


Gas mass [Mq] 


2.9x10'' 


(c) 


4x10'' (d) 


Luminosity [Lq] 


1.9x10'° 


(a) 




Total SFR [Mo yr-'] 


0.4-0.8 


(e) 


0.35 (f) 


Max V„, [km s"'] 


260 


(g) 


200 (h) 



(a) Appleton et al. ( 1981i; (b) Hodge & Miller ( [1992] : (c) Karac hentsev 
& Kashibadze (2006 1; (d) van der Kruit ( [1990) : (e) Gordon et al. (l2004t : 
(f) Walterbos & Braun (1994); (g) Rohlfs & Krertschmann ( fT980t : (h) 
MD05. 



galactic metallicity: M33 is metal poor (its metallicity is similar 
to that of the LMC, Leisy & Dennefeld 2006) and M81 is closer 
in metal contents to the Galaxy. Since M8 1 is also closer to the 
Galaxy in stellar mass content (see Table 7 for a direct compari- 
son of the main properties of these two galaxies), it makes sense 
to compare these two galaxies directly, to explore the relations of 
the metallicity gradients to the characteristics of the galaxy disks 
and to their evolution. The observations of PNe tell us that the 
gradient of M8 1 is steeper than that of our Galaxy in old stel- 
lar populations; the trend seems to persist in the young stellar 
population, but this needs further confirmation. 

From a theoretical point of view, there are two main reasons 
that could influence the gradient slopes, (1) the rotational ve- 
locity of the galaxy and (2) the different situation of these two 
galaxies within their environment. 



MoUa & Diaz (2005) noted that metallicity gradients depend 
on the rotational velocity of the galaxy (thus presumably on its 
total mass), and on its morphological type: more massive galax- 
ies tend to evolve faster and have flatter gradients than lower 
mass galaxies; for a given rotation velocity, gradients are steeper 
for late type than for early type galaxies. Both effects would im- 
ply a steeper gradient for the Galaxy than for M81, contrary to 
the observations. Some other factor must be at play. 

It is likely that the situation of M8 1 within its group of galax- 
ies influences its metallicity gradient. The tidal H i features near 
M8 1 are consistent with a large-scale redistribution of gas in this 
galaxy. M8 1 is the only massive galaxy of its group, surrounded 
by several small galaxies. The interaction of M81 with the dwarf 
galaxies of its group probably steepens the gradient of the ma- 
jor member due to the effect of stripping gas from the external 
regions of the major companion, while interaction of galaxies 
with more or less the same mass would redistribute the gas in 
both galaxies, and thus flatten the gradient. 

The chemical evolution models of Valle et al. (2005) con- 
sider the effects of gas stripping due to galaxy encounters on 
the star formation rate and the evolution of the metallicity. They 
found that, for a stripping occurring 1-3 Gyr after the formation 
of the galaxy and removing 97% of the gas, the region affected 
by the gas removal has a SFR almost a factor of 10 lower than in 
the model without stripping and the relative metallicity is then 
reduced by about 40%. The metallicity reduction is not strongly 
dependent on the time and duration of the stripping episode, but 
is quite sensitive to the relative amount of gas removed from the 
region. This effect should be more pronounced in the outer than 
the inner galactic regions, due to the proximity of the interacting 
galaxies, and it might explain why M8 1 has steeper metallicity 
gradients than the Galaxy. Unfortunately, the models by Valle et 
al. (2005) are limited to a single radial region, thus this conclu- 
sion is of qualitative nature. A radially-resolved chemical evolu- 
tion modeling taking into account the effects of stripping would 
be necessary to determine the effect of tidal interaction on the 
gradient, and also to compare with other results. 

5. Conclusions 

Hectospec/MMT spectroscopy of a sizable sample of PN and 
H II regions in the nearby M81 galaxy has proven very ef- 
ficient to find chemistry of the young and old stellar popu- 
lations, and to pin point the radial metallicity gradients. We 
were able to detect the diagnostic lines for plasma and abun- 
dance analysis in 19 PNe and 14 H n regions. Their analysis 
indicates that the galaxy is clearly chemically enriched, with 
< O/H >H ii reg. / < O/H >PN=1.8, from MMT spectra where 
PNe and H ii regions were simultaneously acquired. We also 
found that there is a noticeable PN metallicity gradient in oxy- 
gen, with Alog(O/H)/ARG=-0.055 dex kpc and that neon and 
sulfur gradient slopes are within 15% of the oxygen one. The 
MMT sample of H n regions have limited galactocentric distribu- 
tion, thus they are insufficient probes of the metallicity gradient. 
The gradient slope from the combined MMT and GS87 H ii re- 
gion samples is steeper (-0.093 dex kpc"') than that of the PNe, 
possibly indicating an evolution of the radial metallicity gradi- 
ents with time: older stellar population show shallower gradi- 
ents, thus gradients are steepening with the time since galaxy for- 
mation. These results have been compared to their homologous 
for the Milky Way and M33, where similar (yet less marked) 
gradient steepening is inferred. We plan to increase the size of 
the H II region sample, and to extend their gradient study to other 
atoms such as neon and sulfur in the near future, for a more com- 
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plete comparison with the PNe, with the goal of obtaining firmer 
observational constraints for the chemical evolutionary models. 
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Table 4. Temperatures and Abundances Uncertainties 
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Table 5. Galactocentric distances and elemental abundances 
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